Transients due to instabilities hinder Kardar-Parisi-Zhang scaling: a unified derivation for surface 
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We propose a unified moving boundary problem for surface growth by electrochemical and chemical vapor 
deposition, which is derived from constitutive equations into which stochastic forces are incorporated. We 
compute the coefficients in the interface equation of motion as functions of phenomenological parameters. The 
equation features the Kardar-Parisi-Zhang (KPZ) non-linearity and instabilities which, depending on surface 
kinetics, can hinder the asymptotic KPZ scaling. Our results account for the universality and the experimental 
scarcity of KPZ scaling in the growth processes considered. 

PACS numbers: 68.35.Ct, 64.60.Ht, 81.15.Gh, 81.15.Pq 



The dynamics of rough surfaces [[TJ is a subject of high in- 
terest. This is due both to its implications for processes of 
technological relevance |^| ||], and to the interesting instances 
that it offers of extended systems evolving in the presence of 
fluctuations [Q]. A very successful framework for the study of 
rough interfaces has been the use of stochastic growth equa- 
tions for the interface height. Among these, the one proposed 
by Kardar, Parisi and Zhang (KPZ) [|J has played a promi- 
nent role, since in particular it has enabled connections to be 
made with other physical problems, like directed polymers in 
disordered media or randomly stirred fluids [[lj. On the basis 
of a coarse-grained description of surface growth and sym- 
metry arguments, the KPZ equation was initially expected to 
describe the dynamics of surfaces growing, e.g. at the expense 
of a vapor phase, in the absence of conservation laws, and is 
thus expected to be relevant to such diverse physical growth 
systems as electrochemical deposition (ECD) [^] or chemi- 
cal vapor deposition (CVD) [|j|. The generality of KPZ scal- 
ing would be a consequence of the phenomenon of univer- 
sality observed for the scaling properties of rough surfaces. 
However, to date very few experiments have been reported 
which are compatible with the predictions of the KPZ equa- 
tion 1^, 0, ||]. Moreover, mere symmetry arguments do not 
enable a detailed connection with phenomenological parame- 
ters describing specific experimental systems, while detailed 
derivations of the KPZ equation were achieved for discrete or 
continuous theoretical models [^J, only indirectly related with 
experiments. These facts have led to invoking additional ef- 
fects on the same coarse-grained level, such as specific noise 
statistics, non-local effects, etc. [y, in order to account for the 
difference between the observed and the predicted scaling be- 
haviors of rough surfaces. However, a wide range of scaling 
exponents ensued, there being no theoretical argument that 
could identify the correct exponents for a specific growth ex- 
periment. 



In this Letter we study two of the main techniques em- 
ployed in experiments on non-conserved surface growth, 
namely ECD and CVD. These techniques have actually 
played a preeminent role in the study of pattern formation 
JlO|], but only recently have they been shown to provide ex- 
perimental realizations of rough interfaces in the KPZ univer- 
sality class [^, ^j. The asymptotic behavior is in both systems 
preceded by exceedingly long unstable transients, extending 
in e.g. the experiments in [^] for up to two days deposition 
time. This complex time behavior is thus far unaccounted for 
on general grounds by any coarse-grained continuum model. 
Here we start from the constitutive equations of ECD and 
CVD, into which we allow for stochastic forces, following a 
similar treatment to that employed in studies of solidification 
Jll|], step dynamics Jl2| ] or fluid imbibition jl3|]. We show 
that ECD and CVD can both be described within a unified 
framework, which provides a stronger statement on universal- 
ity in non-conserved growth phenomena than that restricted 
to scaling behavior. Moreover, we compute the coefficients 
appearing in the ensuing stochastic interface equation of mo- 
tion (IEOM) as functions of the phenomenological param- 
eters characterizing the corresponding physical growth pro- 
cess. The IEOM features the expected KPZ nonlinearity, but 
also instabilities which can hinder asymptotic KPZ scaling. 
Specifically, for the case of non-instantaneous growth events 
at the surface, the IEOM is a stochastic generalization [ [L?] . |T4| | 
of the Kuramoto-Sivashinsky equation JT^], for which very 
long transients due to instabilities are known to occur [ p"o] ] 
before scaling behavior can be observed. Thus, there is no 
need to invoke additional effects at a coarse-grained level in 
order to account for the difficulty in observing KPZ scaling, 
but rather they are due to the long unstable transients which 
will quite generically occur. Our approach also accounts for 
features of discrete growth models [ 17|] and our conclusions 
are thus expected to apply rather generally for non-conserved 
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growth systems. 

We first consider growth by CVD. A successful model of 
this type of growth was formulated and developed in JT^]. 
A stagnant diffusion layer of infinite vertical extent is as- 
sumed to exist above the substrate upon which an aggregate 
will grow. Particles of an intermediate species [concentration 
c(x, z, t) = c(r, t), where x is the coordinate along the initial 
one-dimensional substrate and z is the growth direction] dif- 
fuse through the stagnant layer. When they meet the surface 
they react in order to stick to the aggregate, this occurring with 
an efficiency measured by a kinetic mass-transfer coefficient 
kjj [|l9|. Additional curvature driven effects, such as surface 
diffusion and evaporation-condensation, can influence the lo- 
cal growth velocity of the aggregate. Moreover, we will take 
into account local fluctuations in the vapor phase and surface 
diffusion currents, as well as in the deposition events, in or- 
der to account for the experimental relevance of fluctuations 
for the morphology of surfaces grown by CVD ^ . We thus 
propose the following stochastic generalization of the deter- 
ministic model of CVD 

dc 



dt 



DV 2 c-Vq, (1) 

k D (c - c° eq - Tk + x) = {DVc - q) • n, (2) 
v • n = n(DVc - q) • n- B\7 2 s k - OV • p, (3) 
c(x, z — > oo, t) — c a . (4) 

In Eq. (|lj), D is the diffusion constant, and the conserved noise 
q(r, i) represents fluctuations in the concentration associated 
with diffusion through the stagnant layer. Eqs. (^) and @) hold 
at any point on the aggregate surface r = s, the symbol V 2 
denoting the surface Laplacian operator and n the local unit 
normal vector. Eq. (||) is a mixed boundary condition which 
relates the diffusion current arriving at the aggregate from the 
stagnant layer with the material wich actually deposits, via 
the kinetic coefficient fcj> The noise term \ represents fluc- 
tuations in the deposition events [|l2"|], c® q is the equilibrium 
concentration for a flat inteface, k is the local mean curvature, 
and r = ryc^n/ (ksT), with 7 the surface tension — whose 
anisotropy will be neglected, i.e. we will consider an amor- 
phous or polycrystalline aggregate — , and T temperature. In 
Eq. v is the local aggregate velocity, 51 is the atomic vol- 
ume of the depositing species, and the coefficient of the sur- 
face diffusion current [[TJ B = D s v s jil 2 /kgT, with D s the 
surface diffusivity and v s the surface concentration of parti- 
cles. Moreover, p is a noise term associated with the surface 
diffusion current. Finally, c a in (Q) is a constant value held 
fixed at the edge of the stagnant layer. Note that the determin- 
istic model of CVD [|l8[ is recovered by neglecting q, p, and 
X in (|l])-(Q). We consider zero-mean, uncorrelated and white 
noise terms. A local equilibrium hypothesis 111, 12 1 then al- 
lows us to determine their variances to be 



(g < (r,t)g*(r',0> 
(p^s, t)pV,f)> 
<x(M)x(s',0> 



2£>c(r, t)Sij8(r - r')6(t - t'), (5) 
2D s v s 5 lj S(s - s')5(t - t'), (6) 
{2c(s,t)/k D )8(s-s')6(t-t').(7) 



Before studying the interface dynamics predicted by model 
([j])-(f|), let us show that it also describes electrochemical de- 
position, under a proper interpretation of the fields and pa- 
rameters appearing. For simplicity, we assume a rectangular 
thin cell in which the two electrodes are made of the same 
metal, the cell being filled with a dilute solution of a salt of 
this metal. In a growth experiment by ECD [Qj, an electric 
field is applied driving the motion of cations towards the cath- 
ode, whereupon they stick via a reduction reaction, leading to 
the growth of an aggregate. The constitutive equations, ne- 



glecting convection of the electrolyte, [21 1 describe diffusion 



of cations (concentration C) and of anions (concentration A), 
together with Poisson's equation for the electric field accross 
the cell. This highly non-trivial system can be somewhat sim- 
plified under the electroneutrality condition [|2[ z a A = z c C, 
where ez c and —ez a are the cationic and anionic charges, im- 
plying 



d t C = DV 2 C, 



(8) 



where!) = (fj, c D a + fi a D c )/(fj, a +/x c ) is the ambipolar diffu- 
sion constant, /i C)0 being the cationic and anionic mobilities. 
The electric field configuration and the anion dynamics are 
implicit in the definition of D and in the boundary conditions 
(BC), which we need to specify. The simplest BC are that the 
anion flux be zero both at the anode and at the cathode, where 
only cations contribute to the aggregate growth. The current 
density at the cathode surface is then P] 



Jn = -[z c D c F/(l-t c )]VC -n, 



(9) 



where t c = n c /(n a + fi c ) and F is Faraday's constant. More- 
over, charge transport at the cathode is an activated process, 
whose balance is described by the Butler- Volmer equation 
0,0] 



J= Jo 



(1-/3)11* 



'C/Ca 



(10) 



where Jo is the exchange current density in equilibrium, j3 is 
a coefficient between and 1 describing the asymmetry of the 
energy barrier related to the cation reduction reaction, C a is 
the initial cation concentration, R is the gas constant, and r\ is 
the overpotential, from which a surface curvature contribution 
r\ s has been singled out. By defining the concentration field 
c = D C C/[D(1 — t c )] (c a and c Q eq are defined accordingly) 
and performing an expansion of Eq. (|l0|) for a small value of 
r] s = n ECD j ECD R K /(z c Fk B ) [||, Eqs. (f])-(g) provide the 
stochastic generalization of the ECD model (^])-(|o|) incorpo- 
rating surface tension and surface diffusion effects at the ag- 
gregate surface. Note that Eqs. (^|) and (|l0|) together amount to 
a mixed BC [ p4| ] on c of the type of Eq. (g), with a kinetic co- 
efficient kf) CD = J e-P z <= Fr i/ RT D{l - t c )/{z c FC a D c ). For 
instance, the r\ — > — 00 limit of completely efficient reduction 
at the cathode leads to an absorbing boundary condition there 
of the type c = c? eq + Yn. 

We are now in a position to study the generic model ([l])- 
(0) and draw conclusions for the two diverse growth systems 
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considered. We follow a similar approach to that in [[IT], |l2j 
p3[ | and references therein. Namely, we first note that in the 
zero noise limit Eqs. ([j])-(0) support a fiat solution c(r, t) — 
c(z, t), moving at a constant velocity V = kn [^(c a — c° eq ) — 
1]. By the use of the diffusion Green function, we project the 
moving boundary problem onto the aggregate surface. Finally, 
we perform a perturbation expansion and a long wavelength 
analysis in order to derive a stochastic differential equation for 
a local deviation ({x,t) (in the frame moving with velocity 
V) from the flat interface solution. Details will be reported 
elsewhere [^|. The results are conveniently classified by the 
value of the kinetic coefficient hp [ 25 1 . 



Instantaneous surface kinetics (kD —* oo)- — 
Denoting by ( k (i) the k-th Fourier mode of ((x,t), in the case 
of infinitely fast reaction kinetics at the interface (absorbing 
boundary condition) the IEOM reads 



d t ( k (t)=uj(k)c k (t) + ^T k [(wo 2 } 



■/3l hsolh (t). 



(11) 



Here T k [f(xj\ denotes the fc-th Fourier mode of f(x), the 
non-linear term in ( jll| ) having the expected KPZ form, and 
(3 k hsorh (t) is an additive noise term whose correlations depend 
on the dispersion relation uj(k) [EfJ]. Note that, in principle, 
the system (0)-(||) has multiplicative noise. However, simi- 
larly to [|l2"||, to lowest non-linear order in ({x,t) the IEOM 
features only additive noise terms. For Eq. (nTh the dispersion 



relation reads u(k) = V\k\(l - d l D k 2 )[l - d /l D + (do/4- 
B/D)k 2 Y/ 2 + D(d%/2 - B/D)k 4 - 3d Dk 2 /l D , where we 
have defined a capillarity length do = Til and a diffusion 
length Id = D/V. As we see, due to the shape of uj(k), the 
IEOM ( pi] ) is non-local in space. This is a reflection of the dif- 
fusional instabilities present in the system JlO[]. For instance, 
in the absence of surface diffusion currents (B = 0), and as 
long as do <C Id, the dispersion relation of Eq. ( p"T| ) is of the 
Mullins-Sekerka type [0, oj(k) ~ V\k\(l - d l D k 2 ). How- 
ever, the IEOM does have additive noise and a KPZ non-linear 
term with coefficient V/2, both facts as expected on general 
grounds |jl]]. Although the behavior of Eq. ( pi] ) at large scales 
is not completely known, both a scaling argument and prelim- 
inary numerical simulations indicate that KPZ scaling is not 
asymptotic under these growth conditions p7| ] . 

Non-instantaneous surface kinetics (kD < oo). — 
For the case of non-instantaneous deposition events at the sur- 
face, the dispersion relation oj(k) turns out to depend only on 
even powers of k, which allows to express the IEOM directly 
in configuration space, featuring only local terms [we omit the 
r and t dependencies of £(r, t)]: 

d t ( = -a 2 V 2 ( - a 4 V 4 C + ^(VC) 2 + /3 mixcd (r, t), (12) 

where a 2 = fcc^A, a 4 = k D l D do^/[l - (rfo/'n) 1 ^ 2 ] 
+ B(l + k D /V), and A = 1 - d /l D . Again, /3 mixcd (r,i) is 
an additive noise with w(fc)-dependent correlations p6|]. As 
in Eq. (JTTJ), the coefficient of the KPZ-nonlinearity is V/2. 



However, while ([FT]) always has a band of linearly unstable 
modes, this only happens in Eq. ( |l2"| ) if A > 0, i.e., when 
surface tension is unable to counteract the diffusional insta- 
bilities (do < Id)- In this unstable case, ( p"2| ) is the stochas- 
tic Kuramoto-Sivashinsky (KS) equation, already encountered 
in other interface dynamics contexts, such as step dynamics 
p"2]], or surface erosion by ion-beam sputtering [jl3]. In the 
stochastic KS system, there is a linearly most unstable mode 
k m = (2/£>do) -1 / 2 (assuming B = 0), whose onset time is 
lu^ 1 = 2/ (ko^k^). For asymptotic times, KPZ scaling is 
obtained, but only after an exceedingly long transient [|l6|] . 

From the above analysis of Eqs. ( |TTj ) and (|l2[), KPZ scal- 
ing should be expected only for slow surface kinetics, and 
will nevertheless be affected by early time instabilities, un- 
less the capillary length is larger than the typical diffusion 
length in the system. KPZ scaling corresponds to the confor- 
mal growth mode, which is most interesting for applications 
of CVD-grown films and can indeed be achieved under indus- 
trial conditions by tuning the relative values of kD, Id, and do 
[S, [| . In the case of ECD, we can further verify the predic- 
tions from model (|l])-(0) via Eqs. (pi] ), ( p"2| ) with both experi- 
mental [Q, |28|] and Monte-Carlo [|17|] studies. For instance, in 
the ECD experiment in [28 1 the diffusion length Id ~ 2 cm is 
the largest length scale in the system, and a Mullins-Sekerka 
dispersion relation is reported; consequently, after a stable 
transient the diffusional instabilities that occur completely 
override any scaling behavior at long times. As a difference, 
in the experiments of [Q[ the growing aggregate undergoes an 
unstable transient, beyond which its surface stabilizes into the 
KPZ stationary state. This behavior is qualitatively compati- 
ble with that of the noisy KS equation: from the experimental 
value of the branch spacing one has k rn ~ 1.3 x 10 3 cm -1 . 

2 x 1CT 4 cm s _1 , 
" cm, 



Using 



D ~ 10" 



and V 



one obtains a typical value do — 1/ (2lDk 2 n 



5 x 10" 6 



hence indeed do <C Id — 0.05 cm. Moreover, the instability 
occurs after l/u) m ~ 0.3 x 10 3 s, which allows to estimate 
kD — 8 x 10 -8 cm s _1 V, hence conditions are in the 
slow kinetics regime. 

Eqs. (|])-(Q) also apply to the dynamics of discrete growth 
models such as the MBDLA model [17|, which is a general- 
ization of the DLA model [[j], [l(J to the case of a finite concen- 
tration of random walkers performing biased random walks 
(bias parameter p), which stick to the growing aggregate with 
a finite sticking probability, s. MBDLA reproduces ECD ex- 
periments for one-dimensional substrates quantitatively Jl7|], 
and ECD experiments for two-dimensional substrates quali- 
tatively [5^]. Thus, the sticking probability s plays the role 
of a noise reduction parameter, in the sense that for small s 
values the system reaches faster its KPZ asymptotic scaling 
behavior; this role is played in our continuum model by the 
kinetic mass transfer coefficient fc^ [^9|]. Furthemore, the 
bias p is proportional to the aggregate velocity V. As the 
characteristic branch spacing is, according to (p~2|), approxi- 



mately equal to (lode 



\-l/2 



the continuum model predicts a 



ramified-to-compact transition as p increases, such as is ob- 
served in MBDLA [|l7p. Moreover, since the coefficient of 
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the KPZ term also increases, the corresponding scaling is ex- 
pected to occur earlier as V (or p) increases, again as is ob- 
served in the discrete model. 

In summary, our study of CVD and ECD leads us to expect 
diffusional instabilities to generically hinder KPZ scaling in 
non-conserved growth experiments. Our continuum approach 
from constitutive equations allows to perform detailed com- 
parison with phenomenological parameters in experiments, 
and also provides a physical interpretation for features of dis- 
crete models of kinetic roughening such as noise-reduction 
parameter s [|lj ], e.g. of the type of the sticking parameter in 
MBDLA ||T7j^ More detailed predictions from model (|l])-(0) 
would benefit from a computationally more efficient formula- 
tion. Work along this lines is currently underway Jzo|]. 
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